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Abstract. Energy inequalities are derived for an elliptic-hyperbolic operator 
arising in plasma physics. These inequalities imply the existence of distribution 
and weak solutions to various closed boundary-value problems. An existence 
theorem is proven for a related class of Keldysh equations, and the failure of 
expected methods for obtaining uniqueness is discussed. The proofs use ideas 
recently introduced by Lupo, Morawetz, and Payne for a generalized Tricomi 
operator. The existence of strong solutions under open boundary conditions 
is also proven. MSC2000: 35M10, 35D05, 82D10. 
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1. Introduction 

The equation 
(1) (x - y 2 ) u xx + u yy = 

arises in models of wave propagation through a linear dielectric medium ("cold 
plasma") at frequencies lying below the geometrical optics range; for the physical 
context, see [24]. Here u(x,y), (x,y) G M 2 , is a scalar function. A subscripted 
variable denotes partial differentiation in the direction of the variable. 

The significant property of eq. (pj is that it changes from elliptic to hyperbolic 
type along the parabola 

x - y 2 = 0. 

By analogy with the equations of steady flow, which change from elliptic to hyper- 
bolic type at the speed of sound, it has become conventional to call this parabola 
the sonic curve; in the context of the cold plasma model it is also called a resonance 
curve. Except for a point at the origin, eq. |T|) can be mapped into an equation 
having the same general form as the Tricomi equation 



(2) yU xx + Uyy = 0, 

an equation which is somewhat more accessible than (fTJ). However, both the phys- 
ical and mathematical interest of eq. (|TJ) arise from the tangency of the sonic curve 
to the line x = at the origin. This is the point at which plasma heating might 
occur in the physical model, and a point which appears to be singular in numerical 
studies of solutions; see [IT], [22], and [23] . 

A variety of lower-order terms have been affixed to eq. {T]) in the literature; see, 
e.g., [UJ, [22], and [25]. The variants reflect different toy models for the equation 
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satisfied by the field potential; compare equations (2) and (9) of [22]. The variants 
tend to have the general form 



where the constant k lies in a specified interval. Note that if k = 1, then the 
associated differential operator is formally self-adjoint. 

The formulation of boundary-value problems for eq. ^ is of considerable inter- 
est, as the boundary conditions which are physically natural do not appear to be 
mathematically natural. In particular, it is shown in [17] that the closed Dirich- 
let problem, in which the solution is prescribed on the entire boundary, is over- 
determined for C 2 solutions of (|3|) with n = 1/2. However, the physical properties 
of electromagnetic waves in the cold plasma model suggest that closed boundary- 
value problems are natural and should be correctly posed for eq. (|3j). 

In Sees. 2 and 3 we prove that solutions to closed boundary-value problems for 
([3]) do in fact exist. Although we do not expect classical solutions, we show in 
Sec. 2 the existence of distribution solutions to a homogeneous Dirichlet problem, 
under minimal hypotheses on the domain boundary. In fact these solutions are 
somewhat smoother than conventional distribution solutions, as they lie in the 
function space L 2 . An even smoother distribution solution, which lies in a weighted 
Sobolev space, is derived in Sec. 4 under stronger hypotheses on the boundary. We 
prove the existence of conventional weak solutions to a class of closed boundary- 
value problems, under strong hypotheses on the boundary, in Sec. 3. 

In Sec. 4 we also consider a simpler case in which eq. ([IJ is replaced by a gener- 
alized Cinquini-Cibrario equation [3] 

(4) x 2k+1 u xx + u yy + lower order = 0, 

where k is a non-negative integer. (In fact we do not treat the case k = 0, which 
is the case initially studied by Cinquini-Cibrario.) Equation (j4]) has independent 
mathematical interest as one of two classes of normal forms to which linear, second- 
order elliptic-hyperbolic equations can be reduced near a point [2] • In terms of the 
physical model, eq. ^ would correspond roughly to a resonance surface which 
coincides with a flux surface. In this case the plasma behaves like a perpendicular 
stratified medium and energy absorption occurs along the entire surface, a situation 
more amenable to standard physical arguments than the case in which the two 
surfaces are tangent at a single point; see p. 42 of [23] for a discussion. 

Equation ^ is an example of an equation of so-called Keldysh type — that is, 
an equation of the form 

(5) K(x)u 

XX + Uyy + lower order — 0, 

where K(x) is a continuously differentiable function such that K(Q) — and 
xK(x) > for x ^ 0. Equations of this kind, with various lower-order terms, arise 
in transonic fluid dynamics and singular optics (see, e.g., Sec. 3 of g], an d P3])- 
The arguments of Sec. 4 apply, with only notational alterations, to the slightly more 
general case of a type-change operator having the form K{x) = for n G K, 

n > 3. (The proofs require n to be large enough for K to have three continuous 
derivatives.) 



(3) 
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In Sec. 5 we obtain the existence of unique solutions to eq. (|3]), but under open 
boundary conditions, in which the solution is prescribed on a proper subset of the 
boundary. In Sec. 5.1 we briefly discuss the failure of the preceding methods to 
provide uniqueness in an obvious way. Conditions for the existence of a unique, 
strong solution are introduced in Sec. 5.2. The existence of weak L 2 solutions to a 
large class of open boundary-value problems is proven in Sec. 6 by extending the 
arguments of [18] . 

The approach taken in Sees. 2-4 is modeled on recent work by Lupo, Morawetz, 
and Payne [T3], [5U] on the class of equations 



where K(y) is a continuously differentiable function for which K(0) — and 
yK(y) > for y ^ 0. In the special case K(y) — y, <j6j) reduces to the Tricomi 
equation. For this reason, equations of the form ([6]) are said to be of Tricomi type. 

We adapt the ideas of [13] to the cold plasma context by proving inequalities for 
([3]) having the form 



where C is a positive constant, U and V are function spaces, and L is the differential 
operator of ([3]) with adjoint L*. These energy inequalities are used to show the 
existence of a solution to boundary- value problems in an appropriate function space 
(see, e.g., pQ, Ch. 2). However, we state the inequalities as theorems, rather than 
as lemmas, and derive the existence theorems as corollaries. We do this because in 
their exploitation of inequality (0, the arguments for eqs. ([3]) and ([6]) are essentially 
the same. (While differential operators of Tricomi type are formally self-adjoint, the 
extension of the existence arguments in [13) to the non-self-adjoint case is standard; 
see, e.g., the existence proof in [TS] for weak solutions to a Dirichlet problem for 
([3]) under open boundary conditions.) The arguments for the two equations differ, 
however, in their derivations of the energy inequality itself, which depend on the 
form of the type-change function K. One of the main problems of this paper is to 
find multipliers which allow the Friedrichs abc method to be applied in the right 
way, either in its original form or in the more recent integral variant introduced 
by Didenko [5 . Another is to establish a priori restrictions on the domain which 
allow the method of energy inequalities to be applied. (In Sec. 5.2 we adopt a 
third approach, also due to Friedrichs, in order to establish sufficient conditions for 
uniqueness.) 

In addition to its physical interest, the existence of solutions to closed boundary- 
value problems for equations of the form ([3]) and ([4]) has purely mathematical 
interest as an extension of the methods introduced in |13] to equations which are 
not of Tricomi type. The existence of solutions to closed boundary-value problems 
for certain elliptic-hyperbolic systems which are not of Tricomi type is shown in [7 
and [14) under special conditions on the boundary that do not require the methods 
of [T3]- Our results suggest that the existence — but not necessarily uniqueness — 
of distribution and weak solutions for equations of Keldysh type can be shown by 
arguments closely modeled on those of [13] . 

The following hypotheses on the domain fl are assumed throughout: It is a 
bounded, connected domain of IR 2 having piecewise smooth boundary dfl, oriented 



(6) 



K(y)u xx + Uyy = 0, 
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in a counterclockwise direction; the domain includes both an arc of the sonic curve 
and the origin of coordinates in R 2 (so that, for example, eq. ([3|) is elliptic-hyperbolic 
but not equivalent to an equation of Tricomi type). We neither require nor exclude 
the existence of characteristic arcs on the boundary. We use the term elliptic 
boundary to refer to points [x, y) of the domain boundary on which the type- 
change function K (a;, y) is positive. Similarly, by the hyperbolic boundary we mean 
boundary points for which the type-change function is negative. 

2. Inequalities leading to distribution solutions 

The function spaces introduced in [5] and [T2] reappear in this paper with 
replaced by a different weight function, also denoted by K. In particular, we define 
the space H^fl; K) to be the closure of C^°(f2) with respect to the norm 



(\K\ul + ul + u 2 ) dxdy 



-l 1/2 



where \K\ — \x — y 2 \ . We can write the Hq(Q; jY)-norm in the form 

1/2 

IMIffjCni*) =\J J (\K\ul + ul) dxdy 
as a consequence of the weighted Poincare inequality 



(8) 



MIl 2 (q) - c 



f [ (1*1 

J Jo, 



u y J dxdy. 



Here and below we denote by C generic positive constants, the value of which may 
change from line to line. 

The complexity of the existence arguments is not increased if we replace ([3]) by 
the inhomogeneous equation 

(9) Lu = /, 

where / is a given, sufficiently smooth function of (x, y) and 



(10) 



2^ 

^ 9a; 2 dy 2 dx 
By a distribution solution of equations ([9]), (|10|) with the boundary condition 



(11) u(x,y) = 0V(x,y) e dQ 

we mean a function u £ L 2 (il) such that V£ S H^tt; K) for which L*^ G L 2 (fl), 
we have 



(12) = </>0- 

Here ( , ) denotes the L 2 inner product on fl and ( , ) is the duality bracket associated 
to the H~ l norm [8] 



\ w \\H- 1 (fi;K) 



\(w,Q\ 
sup TTTTi ' 
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Such a solution is a little smoother than the usual notion of distribution solution, 
in which the solution fails to lie in a true function space. 

Theorem 1. Every u £ Cq(£1) satisfies the inequality 

\M\ H i(Q;K) < C\\Lu\\ L 2 (n) , 
where L is defined by with n G [0, 2] , and K = x — y 2 . 

Proof. We consider two cases. 

Case 1: 1 < k < 2. Let 8 be a small, positive constant. Define an operator M 
by the identity 

Mu = au + bu x + cu y 
for a = —1, c = 2 (25 — 1) y, and 

f exp (2^/Qi) if (x,y) e tt+ 
\ exp {65K/Q 2 ) if (re, y) G ' 

where 



0+ = G n\K > 0} 

and f2~ = Choose Qi = exp (25fj,i) , where 



/ii = max -ftT. 



Then V (x, y) G fi + , we have 



2<5if < 25/ii < 2^ ie 2,5A11 = QilogQi. 
Dividing by Qi and exponentiating both sides, we conclude that b < Q\ on + . 
Define the negative number /it 2 by 

/i2 = nun K 

(x, y )en- 

and let Q2 = exp (/x 2 ) ■ Then < Q 2 < 1 and, for given f2, we can choose 5 to be 
so small that 65 < Q 2 . In that case, V (x, y) G O - , 

66^ > 6J^ 2 = 65 log Q 2 > Q 2 log Q 2 . 

We conclude that b > Q 2 on f2~. 

We will estimate the quantity (Mu, Lu) from above and below. Integrating by 
parts, we have 



where 



(Mu, Lu) — J J Ti dxdy 



ti = (auKu x ) x - a x Kuu x - i (aw^ + ^u 2 - aKu 2 x ; 
t 2 = (auuy) - a y uu y - au 2 ; 
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T3 = ±(bKul) x -±(b x K + b)ul; 



T 4 = (bU X Uy) y - byU X Uy 

r 5 = {cUyKux) 



1 



"x 2 



r (ciTtt 2 .! 



c x Ku x u y 
(%K-cy) ul; 



1 r 



T6 = 



i cu l) y - c y u2 y 



the lower-order terms arc: 



T7 



= 9 [( a " 2 )x ~ 



t 8 = 

Tg = KCU y U x . 

As in the Tricomi case considered in [13], one of the coefficients in Mu fails to 
be continuously differentiable on all of £1. When integrating this quantity, a cut 
should be introduced along the line K = separating f2 + from f2~. The boundary 
integrals involving a, 6, and c on either side of this line will cancel by continuity. 

We find that the boundary terms vanish by the compact support of u, and 



(Mu,Lu) = 



n+un- 



uju 2 + au x + 2(3u x u y + ju 2 dxdy, 



where to = 0: 



for 



k - 2 ) b ~ c y> 



and 



<*\n+= [2- — )6K + 2(l-26)y z 



*\n- = [3^--2)5\K\ + 2(l-25)y 2 



0=-{c(K-l)-by] 



1 



7=2^ 



y [26 (6/Qi) + (« - 1) [26 - 1)] < |y| in Q+ 
y [66 (b/Q 2 ) + (k - 1) (2<J - 1)] < n\y\ in Q" 

2 (1 - 6) + * (6/Qx) inO+ 
2(l-6) + 35(b/Q 2 ) infT ' 



On for any scalars £ and i], we have by Cauchy's inequality 

2f3Z V <y 2 e+V 2 , 

so 

< 2 + 2/3^77 + 7?y 2 > a? - (y 2 £, 2 + rj 2 ) + jr] 2 = 



2 - 



5K+(l- 4(5) y 2 + « 



(1 - 25) + 



66 
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> S {Ki 2 + rf) 
provided S is sufficiently small. On Q~, 



so 



2P^<2\ K y^\ <2(y 2 e+V 2 



a£, 2 + 2(3^ + 7 7? 2 > a 2 £ 2 - 2 (y 2 ? + rf) + 1V 2 = 



3 A _ 2 ) S\K\ A6y 2 + L \ 1 b 



e + 5 ( 3— - 2 ) v 2 



>5{\K\e + if)- 

In estimating the coefficient of rj 2 : we used the fact that b > Q2 on f2 _ . In estimating 
the coefficient of £ 2 , we used the fact that 

« - 5) & - 4<% 2 > I - 4<5y 2 > ^ - 4<5y 2 , 

which exceeds zero for 5 sufficiently small with respect to exp (— \1i2D and maxjy?/ 2 . 
Integrating over each subdomain, we obtain 



(13) 



(Mu,Lu)>8 / (\K\u 2 x + u 2 y )dxdy = 5\\u\\ 2 K) . 



We want to obtain an upper bound for (Mu, Lu) in terms of the L 2 -norm of Lu. 
First we recall that max^ ^gn b < Q\ and estimate 



(Mu,Lu)< j I \u\\Lu\ dxdy+ 



(14) C / (Q 1 \u x \ + \y\\uy\)\Lu\dxdy = ii + i 2 . 

J Jn 

The Schwarz inequality and inequality {8} imply that 

h < \\u\\L^n)\\Lu\\ L 2 m < C\\u\\ H i {n . K) \\Lu\\ L 2 {n) . 



Similarly, 



*2 < C 



(\u x \ + \u y \) 2 dxdy 



1/2 



L 2 (fi) 



< C'\\ u \\Hi(n)\\Lu\\ L 2 m . 

We claim that we can choose 6 to be so small that 6/2 > 5 on fi. On f2 + this is 
obvious, as b > 1 there. On fl~ it is almost as obvious, as b > Q2 > 6S > 26. Thus 
we have the additional estimate 



a > S\K\ + 2>$, 



implying that 
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(15) (Mu,Lu)>6\\u\\ 2 H , (n) . 

Substituting the estimates for i\ and ii into inequality (fl4|) and combining this 
with (fT5|) , we obtain 

<%ll^(n) <(Mu,Lu)< 

(16) C (\\u\\ H i (n . K) + \\u\\ ff i in) J \\Lu\\ L 2 m < C'\\u\\ H i (n) \\Lu\\ L 2 {n) . 
Dividing eq. (|16p by the i/p-norm of u, we find that 

Nli?i(n) < C||Lu|| L 2 (0) . 
This completes the proof for the case k £ [1,2], as the norm on the left can 
be replaced by the corresponding weighted norm. (Alternatively, substitute this 
estimate into the extreme right-hand side of (fTr?)) and replace the extreme left-hand 
side of (fT()|) by the extreme right-hand side of (fill)) , to obtain 

\M\ 2 Hi (n; K ) <C||L u ||i 2(n) ; 
then take the square root of each side.) 

Case 2: < K < 1. Again subdivide the domain into fi + and fl~ by introducing 
a cut along the curve K = 0. Integrate by parts as in case 1, choosing a = — 1; 



6 = 

where A is a constant satisfying 



-AA in 
NK in fi" 



(17) 1 + < * < A < l ~ 5 

3 — K AC + 1 

for a sufficiently small positive constant 5, and 

c = -4Ay. 

The boundary integrals involving a and c on either side of the curve will cancel and 
the boundary integrals involving b will be zero on the curve. 
On fi+, 

a= [1 - (1 + k)N] A + 4Ay 2 ; 
/3 = AT(l-2Av)t/; 



7=1 + ^A . 



Also, 



as k G [0, 1), so 



2/3^ > -2A \y£ V \ > -A (y 2 e + if) 
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a£ 2 + 2(3^r, + 1V 2 > a£ 2 ~ N (y 2 £ 2 + rf) + 1V 2 = 
{[1 - (1 + «) N] K + 3Ny 2 } Z 2 + fl + v 2 > ~5 (K£ 2 + V 2 ) , 



by condition (|17p . 
On Or, 



Because 



a = - [I - (3 - k) N] \K\ + Wy 2 ; 
(3 = N(3-2 K )y; 



7 = 1 + -N. 



we have 



2^r)>-6N\y^\, 
af + 2/3^77 + 7 ?7 2 > 



{- [1 - (3 - k) N] \K\ + ANy 2 } e - 3N (y 2 ? + V 2 ) + V 2 

{- [1 - (3 - k)N] \K\ + Ny 2 }e +(l-^ f] 2 >~S (|A|£ 2 + V 2 ) 

again using (fl"7|) . 
We conclude that 



(18) {Mu,Lu)>~5\\u\\ 2 Hh(n . K) . 

In order to obtain an upper bound for (Mu, Lu) in terms of the L 2 -norm of Lu, 
we estimate 

(Mu,Lu)< f f {\u\+N[\K\\u x \+A\y\\u y \]}\Lu\dxdy 



u\\Lu\dxdy + N / [\K\\u x \ + A\y\\u y \]\Lu\ dxdy 
n J Jn 

(19) = H+i 2 . 

As in case 1, the Schwarz inequality and ([8]) imply that 

h < \\u\\mn)\\Lu\\ L 2 {n) < C\\u\\ H i in . K) \\Lu\\ L 2 (n) . 
Because fi is bounded, we can fit it inside a rectangle of the form 



(20) R = {Or, y)\ lR <x< 6 R , p R <y< a R ] . 

for sufficiently large values of \a R \, \0 R \, \jr\, and \Sr\. Define 
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T = 



max 



{|^|+a 2 fl ,l}. 



Then 



i 2 <N 




Vt^\K\\u x \ + A\y\\u y \ \Lu\dxdy 



< N 




Tdxdy \\Lu\\ L 2( n ), 



where 



T = T\K\\u x \ 2 + 8s/T\a R \y/\K\\u x \\u y \ + l6a 2 R \u y \ 2 . 



Applying Cauchy's inequality to T, we obtain 

«2 < C\\u\\ H i {n . K) \\Lu\\ L 2 {Qh 



where C depends on N and R. 

Substituting the estimates for i\ and ii into inequality (|19p . we obtain the desired 
upper bound for (Mu, Lu). Combining this with the lower bound (|18p and dividing 
through by the H^(Cl; K) -norm of u, we complete the proof of Theorem 1. 

Corollary 2. The Dirichlet problem |]|), ilO\) . with k £ [0, 2] possesses a 
distribution solution u £ L 2 (fl) for every f G K). 

Proof. The proof for the case k = 1 is identical to the proof of Theorem 2.2 of 
[13] . with Lemma 2.1 of [13] replaced by Theorem 1 of the present communication. 
Briefly, we define for £ S a linear functional 



and by Theorem 1 (applied to the second term on the right). Now standard Hahn- 
Banach arguments extend the functional to one defined on all of L 2 . The Riesz 
Representation Theorem then guarantees the existence of the distribution solution. 

If k 1, then L is not self-adjoint. Because the operator adjoint to L has the 
form 




by the inequality 



(21) 



\(f,0\< 11/11 



(22) 




estimating L for n in [0, 2] will also yield estimates for L*. Applying the preceding 
argument to the adjoint operator completes the proof of Corollary 2. 
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3. Inequalities leading to weak solutions 

Applications to plasma physics of mixed boundary-value problems are discussed 
in [5]. But we adopt mixed boundary conditions in this section mostly in the 
interest of mathematical generality. The main interest for these equations is in 
closed Dirichlet and Neumann problems. The arguments of this section will extend 
the results of [18 , for an open weak Dirichlet problem for eq. ([I]), to a class of 
generalized, closed Neumann problems. 

If u = (ui,U2) and w = (101,102) are measurable vector- valued functions on $7, 
then the inner product for the space L 2 (Jl; ]R 2 ) will be written 

(u,w) i2(n . R2) = / / (u 1 w 1 +u 2 w 2 )dxdy. 



n 



Analogously, we define H.ic to be the Hilbert space of measurable functions on fi 
for which the weighted _L 2 -norm 



m\H,c 



r r l 1/2 

/ / (\K\ul + ul)dxdy 
J Jn 



is finite; this norm is induced by the the weighted L 2 inner product 



(u, w) K = / / (\K\u1W1 + u 2 w 2 ) dxdy. 
J J n 



In the notation for these spaces, K denotes a diagonal matrix having entries \K\ 
and 1. 

By a weak solution of a mixed boundary- value problem in this context we mean 
an element u £ Ttic(fl) such that 

( 23 ) - ( u > £ * w )l2 ( m2) = (f , w) L 2 (n . R2) 

for every function w G C 1 (£!;R 2 ) for which /C _1 £*w G L 2 (fi;E 2 ) and for which 

(24) wi = OV (x, y) e G 
and 

(25) w 2 =0V(x,y) edQ\G, 

where G is a (possibly empty) subset of <9f2. With a view toward providing a first- 
order generalization of eq. we choose the differential operator C to have the 
form 



(26) 



Kd x d y \ , ( k 



d y -3J 10 



Theorem 3. Let G be a subset of Oil and let K — x — y 2 . Define the functions 
b (x, y) — mK + s and c(y) = [iy — t, where /1 is a positive constant, 



m = 



(/1 + 5) /2 in Q+ 
(p - 5) /2 in n- 
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for a small positive number 5, and t is a positive constant such that [iy — t < Vy £ 
f2. Let s be a sufficiently large positive constant. In particular, choose s to be so 
large that the quantities raK + s, 2cy + s, and b 2 + Kc 2 are all positive. Let 

(27) bdy - cdx < 
on G and 

(28) K {bdy - cdx) > 

on dQ\G. Then there exists a positive constant C such that 

(*, C*M^) > C J J (|if|*? + *!) dxdy 

for any sufficiently smooth 2-vector *, provided conditions \21$ , (25\) are satisfied 
on the boundary for w = A4*, where C* is given by (26\) with k = 1 and 

» = ( -Kc t 

Proof. Again the proof closely follows [13] (Lemmas 4.2 and 4.3, and Theorem 
4.4). After integration by parts and an application of the Divergence Theorem, we 
find that 

(*,£*A4) L2(n . R2) = J j (a*? + 2/3* 1 * 2 + 7 *|) dxdy+ 



(29) 



J (^Kc^l - fc*i* 2 - ^l^j dx + y~ Kb ^i + ^c*i* 2 - dy- 

Because b is not continuously differentiable on f2, it is again necessary to introduce 
a cut along the resonance curve x = y 2 . Evaluating the line integral in (f2"9"|) for our 
choices of b and c, using the fact that c is continuous and that the discontinuous 
term in b vanishes on the resonance curve, we find that the boundary integrals sum 
to zero along the cut. Applying the boundary conditions, we obtain 

J {^Kc^l - 6*1*2 - ^*2^) dx + QiH>*? + ifc*i* 2 - 7^*2) dy 

= \ I ^-*? (b 2 + Kc 2 ) (bdy - cdx) 
2 Jg c 

(30) +\ ( T?*i (& 2 + Kc 2 ) K (bdy - cdx) . 

2 Jdn\G b 

The hypotheses insure that the line integrals in Q30p are finite and nonnegative. 
We have 



a= 2 



Kb x - K x b - (Kc) v + b = (5/2) \ x - y 2 \ + y(w-t) + a/2; 
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Because V e > 



it follows that 



7=2 ( c v - b *) = 



- (1/2) (Sy + t) mil+ 
(1/2) (6y-t) inQ- 

(fj, - 5) /4 in fi+ 
+ 6) /4 in fi- ■ 



2/3*i*2 > 



.(lM±l!ir 4 ;_ e *2, 



/I 



(a* 2 + 2/3*i*2 + 7*2) dxdy > 



5 -\K\+y( l iy-t) + ^ 



(\Sy\ + \t\Y 



* 2 



*2 > dxdy 



> C 



f f (\K\9l + * 2 ) dxdy 
J Jn 



for s sufficiently large, e sufficiently small, and S < /1. This completes the proof. 

Define the space TLjc-i(Q) to consist of measurable functions f = (/1, f2) on f2 
for which /C _1 f lies in the space L 2 (f2;R 2 ) , where 



K,- 1 = 



\k\- 1 
1 

Then H^;-i(r2) is a Hilbert space having inner product 



(31) 



(v,w) 



« K -i(!l) 



= (/C- 2 v, 



'L 2 (fi;R 2 ) 



Corollary 4. Under the hypotheses of Theorem 3, there exists for every f suc/i 
that K,~ 1 M T i S L 2 (f2) a liieafc solution to the mixed boundary-value problem S23]) - 
I125\) with C given by eq. H26\) with k = 0. (The superscripted T denotes matrix 
transpose.) 

Proof. Apply the proof of Theorem 3 in [TB] (a/. |13| . Lemmas 4.2 and 4.3, and 
[2"T]). taking into account that in our case, C 7^ C* . Use Theorem 3 of this paper 
to estimate the formal adjoint of £, which is obtained by taking k = 1 in 
Because A4*f? = w, we have 



(\K\^\ + ^l)dxdy= 11*11^ 



(Q) 



I |-A^f ^Hh^o) = H04 1 w| 1x2(^.^.2) 



for 



/C = 





1 



and 
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a a-X 1 f bw\ — CW 2 

M w — ' 



b 2 + Kc 2 \ cKwi + bw 2 

These norms are finite because of the hypotheses on w and on b 2 + Kc 2 . Using 
Theorem 3, we obtain 

lllfcM-Hli^j < (*,£*.M*) L2(0;R2) = (M-yrwjy^, 

= (/CA1 _1 w,/C _1 /:*w) L2(0 . R2) < \\JCM^ 1 w\\ L 2 {n . M 2 ) \\IC^ 1 £*w\\ L 2 in . M 2 ) , 

in which the right-hand side is finite by our definition of w. Dividing both sides of 
this inequality by the quantity ||/CA / l _1 w|| i 2(Q. R 2), we find that 

2 

| I/CM^WI | L 2 (n . K 2) < - I |£ _1 £*w| |L2 (n . R 2). 

Define the linear functional 

J/(£*w) = (f,w) L2(n . R2) . 

Then 

\J f (£*w) | = | (K^MtXM- 1 ™) | L 2 (0;R2) < C||/C- 1 £*w|| L2(n;R 2 ) . 

We conclude that Jf is bounded on the subspace of H-k;- 1 (^) consisting of elements 
having the form £*~w. Extending the operator to a bounded linear functional on 
the entire space, the Riesz Representation Theorem guarantees the existence of an 
element v S 7i^-i(ri) for which 

(^JyfW) = ( V ^* w )-H /C _ 1 (n) ■ 
The proof is completed by taking u = — /C~ 2 v and applying (|3T|) . 



Remark. If the vector u is sufficiently differentiable, then we can replace ui 
by u x and u 2 by u v for a scalar function u(x,y) . Formally, we then obtain from 
Corollary 4 the existence of a weak solution to eq. ([1]), under mixed boundary 
conditions consisting of Dirichlct conditions 

(32) u x dx + Uydy = 
on G and co-normal conditions 

(33) Ku x dy — u y dx = 
on the complement of G. 

However, the existence of a suitable domain f2 remains to be demonstrated. 

It is apparent that the conditions on are non- vacuous if we consider one of the 
simplest kinds of piecewise smooth domains, a box with a vertex at the origin. Let 
yo be a positive number, let xo = y$, and let fl be the rectangle formed by the line 
segments 
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I : < x < x , y = y ; 

II : x = 0, < y < y ; 
III : < x < x , y = 0; 
IV : x = x , < y < y . 

On the line segments I and II, both (|27|) and (|28| are satisfied. On the line 
segments /// and IV, (|28[l is satisfied. So we can take G to be a nonempty subset 
of lull. Or we can take G to be the empty set, in which case Theorem 3 guarantees 
a weak solution to the conormal problem on Q,. This problem is hard to solve for 
equations of Tricomi type; see [21] and the comments in Sec. 4 of [15] . In the cold 
plasma case, it is the weak Dirichlet problem which is hard to solve, as it is not 
obvious that there is a domain on which G can be taken to be the entire boundary. 

While the restriction of the boundary arcs to vertical and horizontal lines ob- 
viously simplifies the analysis, it is also clear that the existence argument for the 
mixed or co-normal problem extends to domains more general than a rectangle. 
On the elliptic region, all that is required in order for ([28]) to hold with our choices 
of b and c is that dx/ds and dy/ds both be nonnegative, where ds is the element 
of arc length on the boundary. On the hyperbolic region, all that is required for 
either (|2T|) or (|2"5f to hold is that dx/ds and dy/ds both be non-positive. Thus 
a suitable domain might have the form of a lens about that segment of the sonic 
curve x — y 2 which lies in the first quadrant of the xy-plane. It is required that 
the upper boundary of the lens remain in the hyperbolic region of the equations 
without violating the condition dy/ds < 0. For example, let the hyperbolic bound- 
ary be given by the curve y = x q over the interval x € [0, 1] , for q £ (0, 1/2) ; let 
the elliptic boundary be given by the curve y = x r over the interval x S [0,1], for 
r > 1; and let G be a subset of the hyperbolic boundary. 

4. A SOMEWHAT SMOOTHER CLASS OF DISTRIBUTION SOLUTIONS 

It is also possible to strengthen the result of Sec. 2 in a different direction, by 
placing hypotheses on Q which imply more smoothness on the part of the distribu- 
tion solution than mere square-integrability. 

Following Sec. 2 of |12j we consider a one-parameter family (x, y) of inhomo- 
geneous dilations given by 

i, x (x,y) = (A- Q x,A^/), 
where a,(3,\ G M. + . These determine an associated family of operators 

^xu — u o tpx =u\. 

Denote by D the vector field 

= -axd x - (3ydy. 

\\=i 

An open set S] C K 2 is said to be star-shaped with respect to the flow of D if 
V (xq, yo) G ^ and each t £ [0, oo] we have F t (xq, yo) C fl, where 



(34) 



Du 



Tx ux 



F t {x 0l y ) = {x(t),y{t)) = {x Q e- a \y e- pt ) 
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If k is a given function on fl, we define the space L 2 (17; \ k\) and its dual, con- 
sisting, respectively, of functions u for which the norms 



||«||£ 2 (n;|fc|) 

and 

H u llL 2 (n ; |fc|- 1 ) 



v 1/2 

\k\u 2 dxdy I 



/ |fc| 1 u 2 dxdy 
Jn 



1/2 



are finite; see [13], Sec. 3, for details. 

Denote by v a C 1 solution of the Cauchy problem 

(35) Hv = u in Q 

with v vanishing on c?il\{0, 0}, 



(36) lim v (x, y) = 0, 

and 



(37) Hv = av + bv x + cv y . 

We assume in the sequel that v exists. This assumption is justified whenever the 
following conditions are met: f2 is star-shaped with respect to the flow of the vector 
field V = — (b, c) ; b = mx and c = fiy; fi is a positive constant; to is a step function 
in x, taking only positive values, with a single jump at the point x = 0; a is a 
negative constant having sufficiently large magnitude. A proof of the sufficiency of 
these assumptions can be found in step 1 of the proof of Lemma 3.3, [13j . 

Theorem 5. Suppose that x > on Q and that Q is star-shaped with respect 
to the flow of the vector field V = — (b, c) for b = mx and c = fiy, where [i is 
a positive constant and to exceeds 3/i. Then for every u £ C§° ($1) there exists a 
positive constant C for which 

( 38 ) IMlL 2 (n ; |fc|) <C\\Lu\\ H -i {n . k) , 
where k{y) — y 2 and L satisfies ilO\) with k = 1. 

Proof. Let v satisfy eqs. (I3"5l) - ([3"T|) on f2 for a = — M, where M is a positive 
number satisfying 

M= m-3M_ 
2 

for some sufficiently small positive number S. We have the integral identities 

(39) (Iu,Lu) = (v,Lu) = {v,LHv). 

Our choice of the coefficients a, 6, and c are such that a x — a y = b y = c x — and 
all second derivatives also vanish. Substitute into Proposition 12 of the Appendix 
the quantities K = x — y 2 7 n\ = 1, K2 = 0, b = b(x), c — c(y). We have uj = 0, 
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2 7 2 2 

( y - M - <$) x + Sy 2 

P = 0, 



and 



c„ b„ a — m 

1 = -a - Y + Y = M - — =m-2fi-S. 

The boundary integral of Proposition 12 vanishes by the compact support of u. We 
find that if 5 is sufficiently small relative to m and /z, then 

(v,LHv)>5 [ [ (y 2 v 2 x + v 2 y )dxdy. 
J Jo 

The upper estimate is immediate. One applies inequality (|2ip to obtain 

(v,Lu) < \\v\\ H i {n . k) ||iu|| H - l(n .fc) . 

from which the desired inequality follows by the continuity of H as a map from 
Hq (fi; k) into L 2 (f2; |fc|) . This completes the proof of Theorem 5. 

Corollary 6. Let Q be star-shaped with respect to the flow of the vector field 
— V = (mx,fj,y) , where m and fi are defined as in Theorem 5. Suppose that x is 
nonnegative on and that the origin of coordinates lies on <9f2. Then for every 
f G L 2 (yi; l^l" 1 ) there is a distribution solution u € Hq (O; k) to the Dirichlet 
problem (0), where k — y 2 and « = 1. 

Proof. The proof mirrors the arguments for the existence of a distribution solu- 
tion in the proof of Theorem 3.2 in [13j . so we will again be brief. Define a linear 
functional J/ by the formula 

jf{Lt) = u,i). 

Using the fact that L is self-adjoint for k = 1, we estimate 



\Jf (L0\ < C\f\\ L 2 mk \-^\\L(i\\H-Hn-,k)- 
Thus Jf is a bounded linear functional on the subspace of (fi; k) consisting 
of elements having the form L£ with £ 6 (51) . Extending Jf to the entire 
space, the Riesz Representation Theorem guarantees the existence of an element 
u € Hq (57; k) for which 

where £ € -ffg (fo;k) . This completes the proof. 

Remarks, i) While the solution guaranteed by Corollary 2 is only in L 2 , the 
solution guaranteed by Corollary 6 has a derivative in a weighted L 2 -space. Thus 
the solution of Corollary 6 is closer to a conventional weak solution than is the 
solution of Corollary 2. 
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ii) The estimates used to obtain inequality ()38j) extend in a formal way to the 
weight function x — y 2 if, in eq. (|10p . we take k < with |k| sufficiently large. 
Choose 

2k + <5k, 



1 — K \ 2 



4<5 

o = x + mJ(, 

1 — K 



and 



for 



rn 



28 
1 — «* 



5 in 0+ 
—<5 in f2 _ 



K = x — y 2 , and <5 > 0. In that case, w = 0; 



a 



[8{1 + 2\k\)]K + 4Sy 2 in Q+ 
-8K + A5y 2 in Q- ' 



1 r/i ^ , t f 2<5y in ft+ 



and 



The cross terms satisfy 

mn > -2<5 ( y 2 e + v 2 ) , 

so 

/ / (ae + 2^ V + 1V 2 )dxd y >C [ [ (\K\e + V 2 )dxdy 
J Jn J Jn 

provided |k| is sufficiently large. On the boundary between f2 + and Q~, K = 0, so 
along the cut 6 is a smooth function of x only. If f2 is star-shaped with respect to 
the vector field V = — (6, c) , then 

(&,c)-n = 0, 

where n = (—dy,dx) is the unit normal to 0fi . Thus we can obtain the integral 
identity of Proposition 12 in this case as well. It is not obvious that the extension 
is more than formal, as it is not obvious that the vector field V produces a smooth 
solution of the system (|35l) - (|371) . 
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4.1. An equation of Keldysh type. As we remarked in Sec. 1, if we allow the 
resonance curve to be tangent to a flux surface along an entire interval rather than 
at an isolated point, eq. |T]) can be replaced by an equation of the form (01. We 
add lower-order terms to this equation to obtain 

(40) Lu = x 2k+1 u xx + u yy + cix 2k u x + c 2 u = 0, 

where we take k G Z + and require that the constants c\ and c 2 satisfy c± < k + 1 
and c 2 < with |c 2 | sufficiently large. Hypotheses on the magnitudes of lower-order 
terms are quite common for elliptic-hyperbolic equations, such as (|40|) . which are 
not of real principle type; see the remarks in [19j and compare the conditions on 
the lower-order terms in Theorem 4 of [7] . 

Theorem 7. Let a portion of the line x = lie in Q and let the point (0, 0) lie 
on dfl. Assume that ft is star-shaped with respect to the vector field V = — (6, c) , 
where b = rax, and c = \xy. Let \jl be a positive constant and let 

-a/i + n/2£-6/£ in fi+ 
-a/£ + fi/2£ + S/£ in Sl~ 

for a positive constant S, where £ = k + 1 — c%. Let a be a negative constant of 
sufficiently large magnitude. In particular, let a have sufficiently large magnitude 
that to is positive. Then for every w € (Q) there exists a positive constant C 
for which 

\\w\\Li(n;\K\) < C\\L*w\\ H -i {n . K) , 

where K — x 2k+1 and 

L*w = x 2k+1 w xx + w yy + 

(41) (4fc + 2 - ci) x 2k w x + [2k (2fc + 1 - cx) x 2 ^ 1 + c 2 ] w 

is the formal adjoint of the differential operator L of eq. acting on w. 

Proof. The proof is only different in its details from that of Theorem 5. In 
Proposition 12 of the Appendix, take K — x 2k+1 , Ki = (4fc + 2 — c\) x 2k and 
K2 = 2k (2k + 1 — ci) x 2 ^ 1 +C2- Initially we perform all operations over ft + and fl~ 
individually. On the interior of these sub-domains the coefficients are all smooth. 
We find that 

lu = Mx 2k ~ x + c 2 (a - to/2 - fi/2) , 

where M is a constant that depends on k, a, to, and c\ but not on ci- Be- 
cause a is negative and both to and [i are positive, to is positive provided c 2 is 
a negative number having sufficiently large magnitude relative to the quantity 
Mj \a — to/2 — /i/2| . In addition, Proposition 12 implies that 

3 1 

a = -K (a + 2b x ) + - (K x b + Kb x ) + -Kc y - mb 



x 2k+l 



^ - a - (k + 1 - ci) to = S\x\x 2k ; 



20 



THOMAS H. OTWAY 



and 



= \ (K x - ki) c = i ( Cl - 2k - 1) M y:r 2fe ; 



-(l + l/2^)a+(l/4l-l/2)/i-5/2^ in f2+ 
-(l + l/24a+(l/4l-l/2)/i + 5/2^ in O" ' 
Then 7 > £ by the hypotheses on the sign and magnitude of a. We can write (3 in 
the form 

2(3 = (ci - 2fc - 1) ^ya;^ 1 / 2 ^- 1 / 2 . 

Then Ve > 0, 

>-\ci~2k- 1| M |y| fe|x|x 2fc e 2 + -^p^" 1 ^ 2 



£ 



Choose e so small that 



and \a\ so large that 



|ci — 2fe — 1| /imax |y|e < S 



i max |a;| 2fc ~ 1 < - ( 1 + — | a + ( — - - ) a - —. 

en'' V 21) \M 2)^ 21 

We can do this, as |0| is bounded and k > 1. Then, arguing as in the proofs of 
Theorems 1 and 3, we find that 

J J (at 2 + 2(3£rj + ^v 2 ) dxdy >C J J (\K\£ 2 + n 2 ) dxdy. 

The function b(x) fails to be differentiable on the boundary between Cl + and f2~. 
But the coefficients of the boundary terms involving b x either vanish on the line 
K = or cancel out, and the remaining boundary terms are smooth. Thus we can 
integrate over Q + and fl~ , using the support of v and the Divergence Theorem, 
and obtain 



(v,LHv)>C f f (\K\v 2 x + v 2 y )dxdy. 
J Jn 



The remainder of the proof is the same as that of Theorem 5. In particular, the 
existence of a suitable function v still follows from Lemma 3.3 of [T3]: the roles of 
x and y in that argument are symmetric in the sense that the proof is not affected 
if the gap in differentiability is shifted from c(y) to b(x). 

Corollary 8. Let CI and V be defined as in Theorem 7. Then for every f G 
1? (fi; IK]- 1 ) there is a distribution solution u to the Dirichlet problem \4(ty , ill]) 
lying in (ft; K) for K = x 2k+1 . 



Proof. The proof is obtained by arguing exactly as in the proof of Corollary 6, 
substituting the different weight class. 
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Remark. The weight class of Theorem 7 is more natural than that of Theorem 
5, as we expect any singularities that occur to be localized on the sonic curve K. 
However, we do not obtain uniqueness in any obvious way from either theorem, 
because we cannot express weak solutions to either eq. Q or eq. (|40| as the limit 
of an appropriately defined sequence of approximations, as in Sec. 3 of [13) . We 
discuss this problem further in Sec. 5.1. 

5. Inequalities leading to unique solutions 

In the case of equations of Tricomi type, one can say much more. In fact, 
Lupo, Morawetz, and Payne are able to show the existence, with a certain degree 
of regularity, of unique weak solutions to the Dirichlet problem for equations of 
Tricomi type (P2], Sees. 3 and 5). One problem with adapting that approach 
to the present case is the geometry of the sonic curve. The simplest approach is 
to introduce a coordinate transformation which straightens out the sonic curve. 
One then obtains a problem of Tricomi type, but the coordinate transformation is 
singular at the origin, which is precisely the point of interest in the problem. In 
fact, research suggests that the singularity at the origin is an intrinsic property of 
the model [55], [T7], which is not the case, for example, for the Tricomi equation. 
This leads us to expect that solutions to boundary-value problems for eq. Jl} will 
live in a rougher space than solutions to corresponding problems for eq. ([2]). In 
this section we discuss obstructions to obtaining a theorem on the unique existence 
of solutions to closed boundary-value problems for equations of Keldysh type. We 
then impose rather harsh conditions in order to obtain the existence of a unique 
solution to a Dirichlet problem for the cold plasma model. 

5.1. On the difficulty of obtaining uniqueness. A simple but generic example 
will illustrate the difficulty of extending, to equations of Keldysh type, methods 
developed for proving the weak existence of a unique solution to closed boundary- 
value problems for equations of Tricomi type. In our example the energy estimates 
used to obtain Hq (fi; i^)-existence fail for a self-adjoint operator, whereas the 
convergence arguments used to obtain uniqueness seem to require a self-adjoint 
operator. Obviously, our example says nothing about whether a unique solution 
exists, even for this example, but only about the failure of a direct application of 
methods developed for equations of Tricomi type. 
Consider an equation of the form 

(42) [K(x)u x ] x + u yy = 0, 

where the smooth function K(x) changes type on the line x = 0; purely notational 
alterations — e.g., replacing x by (x — xq) in the definition of b(x) — will extend 
the argument to change of type on any vertical "flux line" x — xq. 

In order to guarantee the existence of a smooth solution to the system ([35|) - (|37l) . 
we choose b = mx and c = \iy in that system, where /i is positive constant and m 
is a step function which jumps at the line x = and has the same sign as \i. We 
take a to be a constant of sufficiently large magnitude having the sign opposite to 
that of /i and m. The difficulties will not arise from the choices of \i and a, as the 
argument would be the same for any constants. Moreover, the same difficulties that 
appear for our choice of vector field will appear if b = mx is replaced by b = mf(x), 
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where / is any analytic function vanishing on the flux line (on which m has a jump 
discontinuity). 

The definition of weak solution introduced in [T3] for closed elliptic-hyperbolic 
boundary- value problems is intermediate between the distribution solutions studied 
in the preceding sections and the strong solutions which we will study in Sec. 5.2. 
In accordance with that definition, a weak solution to the Dirichlet problem (|42jl . 
(fTTjl will be a function u £ Hq (O; K) for which 

(Lu,0 = </,0 

for every £ £ Hq (fl;K) , where L is the differential operator of (|42|) . Integrating 
by parts, 



(Lu,t) 



([K(x)i 



'yy 



I ^dxdy 



(Ku x £ x + u y ^y) dxdy. 



In J Jn 

In this case the existence of a weak solution is equivalent to the existence of a 
sequence u n £ C^°(f2) such that 



IK - U\\HZ(SI;K) ~> and \\ Lu n - f\ liT" 1 (fl;K) ~> 

as n tends to infinity. We can thus obtain uniqueness from weak existence by assum- 
ing the existence of two solutions and subtracting their approximating sequences. 

Arguing as in the proof of Theorem 7 in order to establish the existence of a 
solution, we estimate the coefficients w, j3, a, and 7. We find that uj = (3 = 0, and 
that 



K(x) 



Notice that 



a 
2 



xK'(x) 
K{x) 



z r xK'{x) xK'(x) 

(43) hm j—r— = hm - — — — ' ^ = 1. 

v ; x^o K(x) x^o [K(x) - K(0)} 

as K(0) = 0. 

In order for eq. (|42p to change type, K(x) must be monotonic in at least a small 
interval about x — 0. Initially, suppose that K'(x) is positive near x = 0. The 
sum —a + fi/2 cannot be zero, as a and \i have been given opposite sign in order to 
insure the existence of a solution to (|35|) -([37 | . If —a + fi/2 is positive, then for small 
negative values of x, K(x) will be negative. The contribution of [(xK'/K) — 1] m 
will be small by ([4^]) . so a will be negative. If —a + n/2 is negative, then a will be 
negative for small positive values of x for the same reason. An analogous argument 
pertains to the case in which K'(x) is negative near x = 0. 

This example suggests that proving the uniqueness of weak solutions to equations 
of Kcldysh type — and to more complicated non-Tricomi equations — along the lines 
of [13] will require, at the very least, a quite different choice of vector field (6, c), 
which will significantly affect the method. 

For these reasons we will derive uniqueness from the existence of strong solutions 
rather than from the linearity of the differential operator on weak solutions. In 
order to do this, we change the boundary conditions from closed conditions to open 
conditions. 
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5.2. Strong solutions. In the sequel we consider a generalization of the cold 
plasma model: 

(44) Lu = f 
for an unknown vector 

u = (tti (x,y),u 2 (x,y)) 

and a given vector 

f = (fi(x,y) ,h(x,y)) , 

where (x, y) € f2 C M 2 . Here 

(45) (Lu) 1 = [x - a(y)] u lx + Kxm + n 2 u 2yi 

(46) (Lu) 2 = ui y - u 2xi 

where K\ and k 2 are constants; cr(y) > is a continuously differentiable function of 
its argument satisfying 

(47) ff (0) = a'(0) = 0, 

(48) a'{y) > OVy > 0, 
and 

(49) cr'(y) < OVy < 0. 

In the special case in which u(y) = y 2 , k 2 — 0, (/i, /2) = (/, 0) , the components 
of the vector u are continuously differentiable, and u\ = u x , u 2 — u y for some 
twice-differentiable function u(x,y), the first-order system (|44|) - (|46)) reduces to eq. 
©. 

We say that a vector u = (ui,u 2 ) G L 2 (fl) is a strong solution of an operator 
equation of the form (|44ll . with given boundary conditions, if there exists a sequence 
of continuously differentiable vectors, satisfying the boundary conditions, for 
which u v converges to u in L 2 and Lu" converges to f in L 2 . 

A sufficient condition for a vector to be a strong solution was formulated by 
Friedrichs [6] (see also [10]). An operator L associated to an equation of the form 

(50) Lu = A l u x + A 2 u y + Bu, 

where A , A 2 , and B are matrices, is said to be symmetric-positive if the matrices 
A 1 and A 2 are symmetric and the matrix 

Q = 2B* - A l x - A 2 y 
is positive-definite, where B* is the symmetrization of the matrix B : 

B* = i (B + B T ) . 
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If L is not symmetric-positive, then we may consider the equation 

(51) ELu = Ei 

for a non-singular matrix E chosen so that EL is symmetric-positive. 
Define the matrix 

(52) = niA\ an +n 2 A* da , 

where n = (ni,ri2) is the outward-pointing normal vector on dfl. Let M{x,y), 
(x, y) € dfi, be a linear subspace of the vector space V, where u : f2 U dft — ► V. 
Suppose that Af{x,y) depends smoothly on x and y. The boundary condition that 
u lie in J\f is admissible if A/" is a maximal subspace of V and if the quadratic form 
(u, (3u) is non-negative on dfl. 

It is sufficient for admissibility that there exist a decomposition 

a = /?++/?_, 

for which the following three conditions hold: 

i) The direct sum of the null spaces for j3+ and (3- spans the restriction of V to 
the boundary; 

ii) the intersection of the ranges of f3+ and /3_ have only the vector u = in 
common; 

Hi) the matrix [i = /3+ — /3_ satisfies 

If these conditions are satisfied, then the boundary condition 

/3_u = Oon dfl 
is admissible for eq. (|44[) and the boundary condition 

w T pl = on dfl 
is admissible for the adjoint problem 

L*w = g in D.. 

Moreover, both problems can be shown to possess unique, strong solutions. 
Write the system (14"4")) -(14T) ]) in the matrix form 

-°)CI + (-)C)» 
™ + ( o 2 )(;;)■ 

We will show the existence of strong solutions to a subclass of equations for the 
operator L. 

Theorem 9. Assume that on the elliptic boundary of £1 



(54) 



bri\ + cn2 > 0, 
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where n\ and ni are components of the outward-pointing normal vector at each 
point of dVL, and where b and c satisfy, for K = x — o~(y), the inequalities 

(55) 2bn 1 - b x K -b + c y K ~ ca\y) > in fi; 
(2&Ki - b x K - b + c y K - ca'(y)) (2cn 2 + b x - c y ) 

(56) — (bn 2 + CK\ — c x K — c — by) 2 > in f2; 

(57) 2c«2 + b x — Cy > in fl; 

(58) bni - cn 2 < on (dQ)~; 
and 

(59) cKn x + bn 2 > Oon (dSl)~, 

where (dQ) is the hyperbolic boundary. Then equation {44]) > with L given by (53\) 
and the Dirichlet condition 

(60) — 1*1712 + "2^1 = 

imposed on the elliptic portion of dfl, has a unique, strong solution on fl for every 
f I 2 . 

Proof. Define the matrix 

E=( b ~f 

Then the operator EL is symmetric-positive by conditions (|55 p -(|57 p . In order to 
show the existence of strong solutions on SI it is convenient to produce a decompo- 
sition of the matrix 



(61) p 



K (bni — cn 2 ) cKn\ + bn 2 
cKrii + bn 2 — (bni ~ cn 2 ) 

On the elliptic boundary, choose 



fa 

and 



Kbni bn 2 
Kent cn 2 



P- = 



-Kcn 2 Kcn\ 
bn 2 —bn\ 

Then = under the boundary condition (|60[> . Moreover, the intersection and 
range of the two matrices satisfy the conditions for admissibility. We have 
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(i* = {bni + cn 2 ) 

so condition (|54| implies that Dirichlet conditions ([60| are admissible on the elliptic 
part of the boundary. 

On the hyperbolic boundary we choose f3 = (3+ and choose /3_ to be the zero 
matrix. Then the matrix fi* is nonnegative on (dQ) by assumptions (|58[) and (|59p . 
Because the other conditions for admissibility are satisfied trivially on (<9fi) , the 
proof of Theorem 9 is complete. 

Remarks, i) In order to show that inequalities ([55]) - ([59]) are non- vacuous, we 
let a(y) = y 2 , «i = k 2 = 0, b = M + NK/2, and c = Ny, where M and N are 
negative constants, \M\ is sufficiently large, and SI C K x M\R~. Then ([55 |) -([57 | 
are satisfied. Moreover, inequalities (|58|) . (l59|) will be satisfied in a canonical basis 
(ill, H2) = {—dy, dx) provided dy/ds and dx/ds are both nonpositive. This suggests 
that, under the canonical choice of basis, the hyperbolic boundary in Theorem 9 
could be a sufficiently thin lens in the first quadrant, the lower boundary of the lens 
lying along the sonic curve. As a particularly simple example, let the hyperbolic 
boundary be the arc of the circle 

(x-lf+y 2 = 1 

connecting the points (0, 0) and (1, 1) . If \M\ is sufficiently large, then ([54]) requires 
only that dy/ds be bounded below away from zero on the elliptic boundary. 

ii) Under the same choice of basis, in the special case u\ = u Xl u 2 — u y , we 
recover condition ([321) from condition (|60|) and condition (|33p from the adjoint 
condition Ku\ri\ + u 2 n 2 = 0. 

Hi) A hidden smoothness assumption is contained in the choice of the component 
f 2 to be zero in eq. (|44p. as that would imply, in the case u\ = u x , u 2 — u y , 
the equivalence of mixed partial derivatives of the solution. Presumably such a 
condition would be violated at the origin, at which point the difference of the 
mixed partial derivatives might be a delta function. If the difference were somewhat 
smoother than a delta function — that is, an L 2 function, then the methods of this 
section could be applied. 

iv) There is a geometric analogy for condition (|54| : Consider a domain which is 
star-shaped with respect to the flow of a given vector field D satisfying (f3~4"l) . Then 
the boundary will be starlike with respect to D in the sense that cmi + (3n 2 > 0, or, 
in terms of the basis used in remarks i) and ii), (3dx — ady > 0, on the boundary 
(c./. [32]). We have avoided imposing the hypothesis that £1 is Z?-star-shaped in 
Theorems 3 and 9, although it would have been possible to do so formally. The 
reason is that equations of the form ((3|) are only interesting if the origin is included 
in the domain, whereas condition b 2 + Kc 2 > of Theorem 3 and ([55]) of Theorem 
9 are problematic if b and c are homogeneous functions passing through the origin. 

6. Weak solutions in L 2 

The fact that solutions to the closed boundary-value problems in Sees. 3 and 4 
lie in spaces in which a weight function vanishes at the origin is a strong restriction 
on their generality. Theorem 9, with the examples given in the remarks following it, 
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demonstrates the existence of strong solutions which lie in L 2 , even at the origin; 
but unfortunately, the boundary conditions in that theorem are open. In this 
section we show that the existence in L 2 of weak solutions to open boundary-value 
problems is easy to obtain for a wide class of boundaries by arguments which are 
similar to those of [18] . 

Define G to be a subset of the non-characteristic portion of the boundary, d£l\T, 
where T denotes the part of the boundary consisting of characteristic lines. Denote 
by VT(fi) the linear space of continuously differentiable functions (wi,w 2 ) on fi, 
satisfying wi = on G, w 2 = on 9fi\{r U G}, 

(62) widx + u>2dy = V(x, y) G T, 
and 

(i*w) 1 = [x - a (y)] wi x + (1 - Ki) w\ + w 2y , 
(L*w) 2 = wi y - w 2x 

in fi. 

We define a weak solution to eqs. (|44j) -(|46 j) with k 2 — 0, under the mixed bound- 
ary conditions 

(63) u\dx + u 2 dy — V(x, y) S G, 

(64) [x-a (y)\ Ul dy - u 2 dx = V(x, y) € dn\{T U G} 
to be any u G L 2 (Vt) such that Vw G W(Cl), 

(w,f) = -(L*w,u) 
under the L 2 inner product ( , ). 

Theorem 10. Let the noncharacteristic boundary of satisfy the differential 
inequality 

(65) ^ > 

dx m + x 

for a sufficiently large positive constant m and a constant t exceeding 1. Take the 
curve G to be the elliptic boundary ofVL. Let the constant k,\ in eq. exceed 1/2 
and let k 2 — 0. Then Vw G W^fi) there exists a positive constant C for which 

(66) l|w|| L2(0) <G||L*w|| L2(0) . 

Proof. Define the functions b = — (m + x) and c = — ty. We will place various 
conditions on m, all of which require that it be sufficiently large in comparison with 
other parameters — «i, t, |0|, |cr| max (n) and |c'| max (o) — as well as with certain 
explicit combinations of these parameters. By the continuity of a, we can choose 
to so large that the matrix 

be 
M = -Ke b 
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is non-singular on CI, where K — x — cr{y). We have 

(i*w,Mw) = J J^Qdxdy + J J Sdxdy, 



where 

In particular, 



Q = aw 1 + 2(3wiW2 + r yw 2 ■ 



a = (cy - b x ) + I ^ - ki I b - \(j'{y)c 



t - 1 



x+ — — o(y) +m ( Ki y 



1\ ty , 



We can choose to so large that a is bounded below away from zero on Ct. Also, 

7 = ( c y - = —3—' 



P = + & y + K i c ) = ~2~> 



and Ve > 0, 



Choose e to be so small that \K\ty\ e < (t — 1) /2 and m so large that to [kj — (1/2)] > 
/e. Then there is a positive constant C for which 

Q dxdy > C J J (wf + u;|) dxdy. 
Applying the Divergence Theorem, we obtain 

S dxdy = 



K) tyw 1 w 2 



to + x 2 \ m + x 



ty 



- (-K) tyw\ + (to + x)wiw 2 + 



Ml 



, 1 vn ~\~ x 9 
{-K) [ tyw x w 2 H — w 1 



m + x 



mi 



- (-K) tyw\ + (m + x)wiw 2 + y w 



dxdy- 
dxdy = 
dy+ 
dx. 



It is not excluded that the hyperbolic boundary may include one or more char- 
acteristic lines r. Repeatedly applying (|62|) to the terms in wiw 2 on the boundary 
integral over T, we obtain 

— / (Kw\ + w 2 ) [tydx — (to + x) dy] = 0, 
2 Jr 

where on the right we have again used (|62[) . and also the characteristic equations 
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c.f. [18], (3.24)-(3.26). 

On G, wi = and the boundary integral reduces to 

w\ [tydx + (to + x)dy] . 
On 9fi\{GU r}, u>2 — and the boundary integral reduces to 

- / w{{-K)[tydx + {m + x)dy). 

1 Jon\{Gur} 

Both integrals are nonnegative by ([65]) . 
We have shown that 

(LV.Mw) > C||w||| 2 

for some positive number C. Because the elements of M are bounded on Q, applying 
the Schwarz inequality to the inner product (L*w, Mw) yields for all w G W and 
some new constant C > inequality (1661) . 

Corollary 11. Under the hypotheses of Theorem 10, for every f € L 2 (fl) there 
exists on a weak solution to the mixed boundary-value problem |^^[ )- |^[ ) ; i63]). 

Proof. Apply the Riesz Representation Theorem as in 15J (c.f. [18]). 

Remarks, i) This class of boundaries suggests the ice-cream cone-shaped Tri- 
comi domains (see, e.g., Sec. 2 of rotated by 90° in the clockwise direction (so 
that the ice-cream cone is lying on its side, with the cone formed by the intersecting 
characteristic lines in the second and third quadrants). This rotation is expected, 
given the similarity of eq. |T]) to the Cinquini-Cibrario equation, in which the sonic 
curve is rotated 90° with respect to the sonic curve for the Tricomi equation J2]). 
In fact, the sonic curve of ([1]) is approximated near the origin by the sonic curve of 
the Cinquini-Cibrario equation. 

ii) More generally, Theorem 10 and Corollary 11 remain true for any choice of 
M for which 

dy > _c 
dx ~ b 

on the characteristic boundary, a and 7 are bounded below by a positive constant, 
and «7 — (3 2 is nonnegative. 

Hi) In the corresponding theorem of [TS], the origin of coordinates was forced 
to lie on the boundary of the domain. The question of whether weak solutions to 
boundary- value problems for eqs. (|44|) - (|46[) can be shown to exist for cases in which 
the origin is allowed to be an interior point was raised in Ch. 3 of |25j . In allowing 
the origin to lie at either a boundary point or an interior point, we have shown the 
answer to that question to be "yes." 
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7. Appendix. A multiplier identity 

The proofs of Theorems 5 and 7 are based on a fundamental multiplier identity 
for the composition of a second-order operator and a first-order operator. The proof 
of the identity is elementary, as it is based on integration by parts followed by an 
application of the Divergence Theorem. But due to its importance, we provide a 
full derivation (in slightly greater generality than we need). 

Define the operator L on functions d€C 3 (O) with v = on dSl, by 

Lv = K (x, y) v xx + v yy + kiv x + k 2 v, 
where the type-change function K and the lower-order coefficients k\ — Ki(x,y), 
«2 = K2 (x, y) are all C 3 functions. Define 

Hv = av + bv x + cv y , 
where a is a constant; c is a linear function of x and y; b is linear in x but possibly 
nonlinear in y; b and c have vanishing mixed partial derivatives. 

Proposition 12. 

/ [ v ■ LHv dxdy =\ £> [Kv x + v 2 ) (cdx — bdy) 
J Jn 2 J gn 

+ J J uov 2 + avl + 2(3v x v y + -fVydxdy, 



where 



2u = (K xx - k 1x + 2k 2 ) a - [(K xx - k 1x + n 2 ) b] a 
- [(K xx - ki x + k 2 ) c] ; 



c 



a= ^ JL _ JL _ a j K+ ^_ Kx _ Kl j b+ _ Ky . 
2(3=(K x -K 1 )c-{c x K + b y ); 

Proof. Writing 

LHv = K [(a + 2b x ) v xx + bv 

+ (a + 2c y ) 

Uyy ~t~ byy^x H~ ^^y^xy ^^xyy H~ C^yyy 

+Ki [(a + b x ) v x + bv xx + c x v y + cv yx ] + k 2 (av + bv x + cv y ) , 

we have 

16 



v ■ LHv — Ti, 



i=i 
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n = vK (a + 2b x ) • 



v xx = | (a + 2b x ) ^ 



(a + 2b x ) ( Kv x - 2 KxV 



-K (a + 2b x ) v 2 x + -K xx (a + 2b x ) v 2 ; 
t 2 = vKbv xxx = 



-^Kbvl + 



b \Kv xx - K x v x + -K xx vj + b x (K x v - Kv x ) 

1 3 

- 2 {Kxxxb + 3K xx b x ) v 2 + - (Kb) x v 2 x ; 



t 3 = 2vKc x v yx = 2 (vKc x v y ) x - (K x c x v 2 ) y - 2Kc x v x v y + K xy c x v 2 ; 

r 4 = vKcv yxx = 

{v [c(Kv yx - K x v v ) - Kc x v y ]} x - ji^f cv 2 x - (^K xx c + K x c x ^j « 2 | 

\ \} Kxx y c + KxxC v + ^K xy c x ] v 2 - (Kc) y ^ j + {Kc) x v x v y ; 
T 5 = v(a + 2c y ) Vyy = [v ( a + 2^) ^ - (a + 2c y ) v 2 y , 
t 6 = vb yy v x = i (b yy v 2 ) x ; 

T 7 = 2vbyV X y = (2vbyV X ) y - 2byV X Vy ~ ^yyV^^, ] 

= VbV X yy = - {^bVy^j + [(bV X y ~ byV X )v] y + ^ [b X + (^J,^) J + byV X 

T 9 = VCVyyy = ~ (ci^ + [(CVyy ~ CyVy) V ] y + ^CyV^, 

no = (a + b x ) v x = 1 [(a + bx) ^i« 2 ] x - (a + b x ) v 2 ; 

Tn = VK!bv xx = 



K\bv x - - {K X b) x v 



v\ + i {n lxx b + 2n lx b x ) v 2 ~ Kibvl; 



[K\C X V j y K\yC x 



T 13 = VK lC V yx = {vKiCV y ) x - - [(K lC ) x V 2 ] y 
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— KiCV x V- 



X<Jy 



KlyC x + (m x c) 



v 2 ; 



T14 = K 2 av ; 



T15 = vk^z = ^ ( K 2&w 2 ) x - ^ ( K 2&) x w 2 ; 

Tig = WK 2 CW a = i (c^'y 2 )^ - ^ (CK 2 )j, W 2 - 

Collect terms and integrate over f2. Applying the Divergence Theorem, taking 
into account that v (but not necessarily v x or v y ) vanishes on <9fi, completes the 
proof. 

Similar estimates are applied, in the proof of Theorem 1, to a product having the 
simpler form (Hv, Lv). The matrix identities which underlie the proofs of Theorems 
3 and 10 are analogous but also simpler, as the differential operators in those cases 
are first-order and the equations are in reduced form. In Proposition 12 the operator 
is effectively third-order and there are zeroth-order terms in the equations. 

Acknowledgment. I am grateful to Profs. D. Lupo and K. R. Payne for dis- 
cussion of this problem, and to the referee for extremely helpful comments. 
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